
	

*Change the directory to the location of the replication files on your computer
cd "C:\PCult_Replication"

do "CreateDataSet"


	summarize TrumpIndex
	*.3920479
	hist TrumpIndex
	tab TrumpIndex TrumpCult

	
	
	 centile Conservatism1 GOP  cOpenness1 cNeuroticism1 cAgreeableness1 cExtraversion1 cConscientiousness1 if TrumpIndex>.3920479 & TrumpIndex~=., centile(5 95)
	 
/*

                                                          Binom. interp.   
    Variable |       Obs  Percentile    Centile        [95% conf. interval]
-------------+-------------------------------------------------------------
Conservati~1 |       475          5    .3333333             .25        .375
             |                   95        .875        .8333334    .9166666
         GOP |       475          5           0               0           0
             |                   95           1               1           1
  cOpenness1 |       475          5    .1796875        .1503906      .21875
             |                   95    .7644531        .7285859    .8000153
cNeurotici~1 |       475          5    .1093986        .0640757    .1390978
             |                   95    .7913533        .7595337    .8363832
cAgreeable~1 |       475          5    .4173685        .3932878    .4332703
             |                   95    .9326315        .9001262    .9607542
cExtravers~1 |       475          5    .2838129        .2293165    .3156475
             |                   95    .7861511        .7377982    .8170934
cConscient~1 |       475          5    .3955517        .3590861    .4092527
             |                   95    .9759787          .95021    .9810071


*/
	. 


	 
	 
**** Figure 2

***Personality Only
  
  set scheme stcolor

*mean	
	logit TrumpCult  cOpenness1 cNeuroticism1  cAgreeableness1 cExtraversion1 cConscientiousness1 if TrumpIndex>.3920479 & TrumpIndex~=.

estat ic
 
margins if cOpenness1>.1796875  & cOpenness1<.7644531  & cNeuroticism1>.1093986 & cNeuroticism1<.7913533 & cAgreeableness1 > .4173685  & cAgreeableness1<.9326315 & cExtraversion1 > .2838129  &  cExtraversion1 < .7861511 & cConscientiousness1 > .3955517 & cConscientiousness1<.9759787, dydx( *) 
 
marginsplot,  title("Study 1 (BFI)") subtitle("") ytitle("") derivlabels  xtitle("") horizontal plotopts(connect(none)) xline(0) xscale(r(-.6 .8)) xlabel(-.6(.4).8) xtick(-.6(.4).8)
 
   
graph save "Graph" "meanS101.gph", replace


***Personality Conservatism1 GOP
  
set scheme stcolor
  
logit TrumpCult Conservatism1 GOP cOpenness1 cNeuroticism1  cAgreeableness1 cExtraversion1 cConscientiousness1 if TrumpIndex>.3920479 & TrumpIndex~=.

estat ic
 
 
 margins if Conservatism1>.3333333 & Conservatism1<.875 & cOpenness1>.1796875  & cOpenness1<.7644531  & cNeuroticism1>.1093986 & cNeuroticism1<.7913533 & cAgreeableness1 > .4173685  & cAgreeableness1<.9326315 & cExtraversion1 > .2838129  &  cExtraversion1 < .7861511 & cConscientiousness1 > .3955517 & cConscientiousness1<.9759787, dydx( *) 
 
marginsplot,  title("Study 1 (BFI and Controls)") subtitle("") ytitle("") derivlabels  xtitle("") horizontal plotopts(connect(none)) xline(0) xscale(r(-.6 .8)) xlabel(-.6(.4).8) xtick(-.6(.4).8)
 
graph save "Graph" "meanS102.gph", replace
  
  
  
***Figure 3 
******Trump INdex Dep Var >0.9
  
gen Trump_Supporter_S1 = 0
replace Trump_Supporter_S1 = 1 if TrumpIndex>.90 & TrumpIndex~=.


tab TrumpIndex Trump_Supporter_S1
summarize Trump_Supporter_S1
  
set scheme stcolor
  
	
logit Trump_Supporter_S1 Conservatism1 GOP cOpenness1 cNeuroticism1  cAgreeableness1 cExtraversion1 cConscientiousness1 if TrumpIndex>.3920479 & TrumpIndex~=.

estat ic
 
 
  margins if Conservatism1>.3333333 & Conservatism1<.875 & cOpenness1>.1796875  & cOpenness1<.7644531  & cNeuroticism1>.1093986 & cNeuroticism1<.7913533 & cAgreeableness1 > .4173685  & cAgreeableness1<.9326315 & cExtraversion1 > .2838129  &  cExtraversion1 < .7861511 & cConscientiousness1 > .3955517 & cConscientiousness1<.9759787, dydx( *) 
 

marginsplot,  title("Study 1 (BFI and Controls)") subtitle("") ytitle("") derivlabels  xtitle("") subtitle("Trump Index > 0.9") horizontal plotopts(connect(none)) xline(0) xscale(r(-.6 .8)) xlabel(-.6(.4).8) xtick(-.6(.4).8)
 
  graph save "Graph" "mean09_S1.gph", replace
  
  
  
***Personality Conservatism1, GOP only
  
set scheme stcolor
  
logit TrumpCult Conservatism1 cOpenness1 cNeuroticism1  cAgreeableness1 cExtraversion1 cConscientiousness1 if GOP==1

estat ic
 
margins if Conservatism1>.375 & Conservatism1<.9166666 & cOpenness1>.1591797 & cOpenness1<.7768555 & cNeuroticism1>.0869362  & cNeuroticism1<.7504699 &  cAgreeableness1>.4184211 & cAgreeableness1<.9585526 & cExtraversion1>.2877698  & cExtraversion1<.807554 & cConscientiousness1 >  .4190392 & cConscientiousness1< .9790925, dydx( *)
 
marginsplot,  title("Study 1 (BFI and Conservative)") subtitle("GOP Only") ytitle("") derivlabels  xtitle("") horizontal plotopts(connect(none))  xline(0) xscale(r(-.6 1.0)) xlabel(-.6(.2)1.0) xtick(-.6(.2)1.0)
 
graph save "Graph" "GOPS102.gph", replace 
 


***Figure 4, Facets 
  
  
  centile cIdeas1_Open  cAesth1_Open cDepress1_Neu cAnxiety1_Neu  cComply1_Agr cAltruism1_Agr cAssertive1_Ext cActivity1_Ext cOrder1_Con cDiscipline1_Con GOP if TrumpIndex>.3920479 & TrumpIndex~=., centile(5 95)

 /*
                                                          Binom. interp.   
    Variable |       Obs  Percentile    Centile        [95% conf. interval]
-------------+-------------------------------------------------------------
cIdeas1_Open |       475          5    .2594228        .2139968    .2918372
             |                   95    .7833616        .7573378    .8302207
cAesth1_Open |       475          5    .2609113        .2150986    .3021584
             |                   95    .8215828        .7891414    .8752998
cDepress1_~u |       475          5    .2305085        .2090396    .2542373
             |                   95    .8531073        .8307117    .8793896
cAnxiety1_~u |       475          5       .0625               0       .0625
             |                   95        .875           .8125        .875
cComply1_Agr |       475          5    .3320099        .2977667    .3696378
             |                   95    .9434242         .930521    .9544713
cAltruism1~r |       475          5    .3446154        .3038462    .3730769
             |                   95    .8853846        .8616768    .9217383
cAsser~1_Ext |       475          5          .2             .15         .25
             |                   95          .8             .75          .8
cActivity1~t |       475          5    .2580645        .2193548    .2965421
             |                   95    .8670968         .851845    .8903226
 cOrder1_Con |       475          5    .3736842        .3421053    .3945475
             |                   95    .8947368        .8842105    .9052632
cDisci~1_Con |       475          5    .3313725        .3099609    .3528824
             |                   95    .9624183        .9462547    .9689543
         GOP |       475          5           0               0           0
             |                   95           1               1           1


*/
  
    logit TrumpCult GOP cIdeas1_Open  cAesth1_Open cDepress1_Neu cAnxiety1_Neu  cComply1_Agr cAltruism1_Agr cAssertive1_Ext cActivity1_Ext cOrder1_Con cDiscipline1_Con if TrumpIndex>.3920479 & TrumpIndex~=.

  set scheme stcolor
  estat ic
 
 margins if cIdeas1_Open>.2594228  & cIdeas1_Open<.7833616 & cAesth1_Open>.2609113 & cAesth1_Open<.8215828  & cDepress1_Neu > .2305085  & cDepress1_Neu < .8531073 & cAnxiety1_Neu >  .0625  & cAnxiety1_Neu< .875 & cComply1_Agr > .3320099  & cComply1_Agr < .9434242   &  cAltruism1_Agr>.3446154 & cAltruism1_Agr<.8853846   & cAssertive1_Ext>.2 & cAssertive1_Ext<.8 & cActivity1_Ext>.2580645  & cActivity1_Ext<.8670968  & cOrder1_Con>.3736842  & cOrder1_Con < .8947368 & cDiscipline1_Con>.3313725 & cDiscipline1_Con < .9624183 , dydx( *) 
 
  
  marginsplot,  title(Trump Followers) subtitle("") ytitle("") derivlabels xtitle("") horizontal plotopts(connect(none))  xline(0) xscale(r(-.4 .6)) xlabel(-.4(.2).6) xtick(-.4(.2).6)
 
 
 
 graph save "Graph" "MeanFacetsTrump.gph", replace
 
  
 *ConservativeTen

 tab ConservativeTen
 
 logit ConservativeTen GOP cIdeas1_Open  cAesth1_Open cDepress1_Neu cAnxiety1_Neu  cComply1_Agr cAltruism1_Agr cAssertive1_Ext cActivity1_Ext cOrder1_Con cDiscipline1_Con if TrumpIndex>.3920479 & TrumpIndex~=.


  estat ic
 
 margins if cIdeas1_Open>.2594228  & cIdeas1_Open<.7833616 & cAesth1_Open>.2609113 & cAesth1_Open<.8215828  & cDepress1_Neu > .2305085  & cDepress1_Neu < .8531073 & cAnxiety1_Neu >  .0625  & cAnxiety1_Neu< .875 & cComply1_Agr > .3320099  & cComply1_Agr < .9434242   &  cAltruism1_Agr>.3446154 & cAltruism1_Agr<.8853846   & cAssertive1_Ext>.2 & cAssertive1_Ext<.8 & cActivity1_Ext>.2580645  & cActivity1_Ext<.8670968  & cOrder1_Con>.3736842  & cOrder1_Con < .8947368 & cDiscipline1_Con>.3313725 & cDiscipline1_Con < .9624183 , dydx( *) 
 
 
  marginsplot,  title(Conservative) subtitle("") ytitle("") derivlabels xtitle("") horizontal plotopts(connect(none))  xline(0) xscale(r(-.4 .6)) xlabel(-.4(.2).6) xtick(-.4(.2).6)
 
 
 
 graph save "Graph" "MeanFacetsCon.gph", replace
 
 
 ***combine
 
  graph combine MeanFacetsTrump.gph MeanFacetsCon.gph, rows(1) cols(2)  title(Facets of the Big Five) name(Facets, replace) saving(Facets, asis replace)
  
 
****Figure 5, kdensity plots
 
 
*cConscientiousness1

set scheme stcolor

kdensity cConscientiousness1  if TrumpCult==1, addplot((kdensity cConscientiousness1  if TrumpCult==0 & TrumpIndex>.3920479 & TrumpIndex~=.)) title("") subtitle("") ytitle("Above Mean", size(medlarge)) xtitle("") legend(off)

graph save "Graph" "kd_ConscMean.gph", replace


kdensity cConscientiousness1  if TrumpCult==1, addplot((kdensity cConscientiousness1  if TrumpCult==0)) title("Conscientiousness") subtitle("") ytitle("All", size(medlarge)) xtitle("") legend(off)

graph save "Graph" "kd_Consc.gph", replace


kdensity cConscientiousness1  if TrumpCult==1 & ConservativeTen==1, addplot((kdensity cConscientiousness1  if TrumpCult==0 & ConservativeTen==1)) title("") subtitle("") ytitle("Conservatives", size(medlarge)) xtitle("") legend(off)

graph save "Graph" "kd_ConscCons.gph", replace


kdensity cConscientiousness1  if TrumpCult==1 & GOP==1, addplot((kdensity cConscientiousness1  if TrumpCult==0 & GOP==1)) title("") subtitle("") ytitle("GOP", size(medlarge)) xtitle("") legend(off)

graph save "Graph" "kd_ConscGOP.gph", replace

kdensity cConscientiousness1  if TrumpCult==1 & GOP==0, addplot((kdensity cConscientiousness1  if TrumpCult==0 & GOP==0)) title("") subtitle("") ytitle("Non-GOP", size(medlarge)) xtitle("") legend(off)

graph save "Graph" "kd_ConscNotGOP.gph", replace



*cOrder1_Con


kdensity cOrder1_Con  if TrumpCult==1, addplot((kdensity cOrder1_Con if TrumpCult==0 & TrumpIndex>.3920479 & TrumpIndex~=.)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_OrdMean.gph", replace


kdensity cOrder1_Con  if TrumpCult==1, addplot((kdensity cOrder1_Con if TrumpCult==0 )) title("Order") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_Ord.gph", replace

kdensity cOrder1_Con  if TrumpCult==1 & ConservativeTen==1, addplot((kdensity cOrder1_Con  if TrumpCult==0 & ConservativeTen==1)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_OrdCons.gph", replace



kdensity cOrder1_Con  if TrumpCult==1 & GOP==1, addplot((kdensity cOrder1_Con if TrumpCult==0 & GOP==1)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_OrdGOP.gph", replace


kdensity cOrder1_Con  if TrumpCult==1 & GOP==0, addplot((kdensity cOrder1_Con if TrumpCult==0 & GOP==0)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_OrdNotGOP.gph", replace




*cDiscipline1_Con



kdensity cDiscipline1_Con if TrumpCult==1, addplot((kdensity cDiscipline1_Con if TrumpCult==0 & TrumpIndex>.3920479 & TrumpIndex~=.)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_DiscMean.gph", replace


kdensity cDiscipline1_Con if TrumpCult==1, addplot((kdensity cDiscipline1_Con if TrumpCult==0)) title("Discipline") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_Disc.gph", replace

kdensity cDiscipline1_Con   if TrumpCult==1 & ConservativeTen==1, addplot((kdensity cDiscipline1_Con   if TrumpCult==0 & ConservativeTen==1)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_DiscCons.gph", replace



kdensity cDiscipline1_Con if TrumpCult==1 & GOP==1, addplot((kdensity cDiscipline1_Con if TrumpCult==0 & GOP==1)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_DiscGOP.gph", replace

kdensity cDiscipline1_Con if TrumpCult==1  & GOP==0, addplot((kdensity cDiscipline1_Con if TrumpCult==0  & GOP==0)) title("") subtitle("") ytitle("") xtitle("") legend(off)

graph save "Graph" "kd_DiscNotGOP.gph", replace



graph combine kd_Consc.gph kd_Ord.gph kd_Disc.gph kd_ConscMean.gph kd_OrdMean.gph  kd_DiscMean.gph kd_ConscCons.gph kd_OrdCons.gph  kd_DiscCons.gph kd_ConscGOP.gph kd_OrdGOP.gph  kd_DiscGOP.gph kd_ConscNotGOP.gph kd_OrdNotGOP.gph kd_DiscNotGOP.gph, rows(5) cols(3)  title("") note(Blue lines: Trump loyalists; Red lines: others in category) name(FigKD, replace) saving(FigKD, asis replace)
  
  
  
 ***Figure 6, Facets, Aspects
 
  
****Trump_JobApproval

logit Trump_JobApproval Conservatism1 GOP cIdeas1_Open  cAesth1_Open cDepress1_Neu cAnxiety1_Neu  cComply1_Agr cAltruism1_Agr cAssertive1_Ext cActivity1_Ext cOrder1_Con cDiscipline1_Con if TrumpIndex>.3920479 & TrumpIndex~=.

  set scheme stcolor
  estat ic
 
 margins if cIdeas1_Open>.2594228  & cIdeas1_Open<.7833616 & cAesth1_Open>.2609113 & cAesth1_Open<.8215828  & cDepress1_Neu > .2305085  & cDepress1_Neu < .8531073 & cAnxiety1_Neu >  .0625  & cAnxiety1_Neu< .875 & cComply1_Agr > .3320099  & cComply1_Agr < .9434242   &  cAltruism1_Agr>.3446154 & cAltruism1_Agr<.8853846   & cAssertive1_Ext>.2 & cAssertive1_Ext<.8 & cActivity1_Ext>.2580645  & cActivity1_Ext<.8670968  & cOrder1_Con>.3736842  & cOrder1_Con < .8947368 & cDiscipline1_Con>.3313725 & cDiscipline1_Con < .9624183 , dydx( *) 
 
  marginsplot,  title(Trump Job Approval) subtitle("") ytitle("") derivlabels xtitle("") horizontal plotopts(connect(none))  xline(0) xscale(r(-.6 .6)) xlabel(-.6(.2).6) xtick(-.6(.2).6)
 
 graph save "Graph" "jobMean.gph", replace
 
 
 ****Trump_BigLie
 
 logit Trump_BigLie Conservatism1 GOP cIdeas1_Open  cAesth1_Open cDepress1_Neu cAnxiety1_Neu  cComply1_Agr cAltruism1_Agr cAssertive1_Ext cActivity1_Ext cOrder1_Con cDiscipline1_Con if TrumpIndex>.3920479 & TrumpIndex~=.

  set scheme stcolor
  estat ic
 
margins if cIdeas1_Open>.2594228  & cIdeas1_Open<.7833616 & cAesth1_Open>.2609113 & cAesth1_Open<.8215828  & cDepress1_Neu > .2305085  & cDepress1_Neu < .8531073 & cAnxiety1_Neu >  .0625  & cAnxiety1_Neu< .875 & cComply1_Agr > .3320099  & cComply1_Agr < .9434242   &  cAltruism1_Agr>.3446154 & cAltruism1_Agr<.8853846   & cAssertive1_Ext>.2 & cAssertive1_Ext<.8 & cActivity1_Ext>.2580645  & cActivity1_Ext<.8670968  & cOrder1_Con>.3736842  & cOrder1_Con < .8947368 & cDiscipline1_Con>.3313725 & cDiscipline1_Con < .9624183 , dydx( *) 
 
  marginsplot,  title(Trump Election Win) subtitle("") ytitle("") derivlabels xtitle("") horizontal plotopts(connect(none))  xline(0) xscale(r(-.6 .6)) xlabel(-.6(.2).6) xtick(-.6(.2).6)
  
graph save "Graph" "BigLieMean.gph", replace
 

 ****Trump_Leader
 
logit Trump_Leader Conservatism1 GOP cIdeas1_Open  cAesth1_Open cDepress1_Neu cAnxiety1_Neu  cComply1_Agr cAltruism1_Agr cAssertive1_Ext cActivity1_Ext cOrder1_Con cDiscipline1_Con if TrumpIndex>.3920479 & TrumpIndex~=.

  set scheme stcolor
  estat ic
 
 margins if cIdeas1_Open>.2594228  & cIdeas1_Open<.7833616 & cAesth1_Open>.2609113 & cAesth1_Open<.8215828  & cDepress1_Neu > .2305085  & cDepress1_Neu < .8531073 & cAnxiety1_Neu >  .0625  & cAnxiety1_Neu< .875 & cComply1_Agr > .3320099  & cComply1_Agr < .9434242   &  cAltruism1_Agr>.3446154 & cAltruism1_Agr<.8853846   & cAssertive1_Ext>.2 & cAssertive1_Ext<.8 & cActivity1_Ext>.2580645  & cActivity1_Ext<.8670968  & cOrder1_Con>.3736842  & cOrder1_Con < .8947368 & cDiscipline1_Con>.3313725 & cDiscipline1_Con < .9624183 , dydx( *) 
 
  marginsplot,  title(Trump Leadership) subtitle("") ytitle("") derivlabels xtitle("") horizontal plotopts(connect(none))  xline(0) xscale(r(-.6 .6)) xlabel(-.6(.2).6) xtick(-.6(.2).6)
  
  graph save "Graph" "leaderMean.gph", replace
 

 ****Trump_Status
 
logit Trump_Status Conservatism1 GOP cIdeas1_Open  cAesth1_Open cDepress1_Neu cAnxiety1_Neu  cComply1_Agr cAltruism1_Agr cAssertive1_Ext cActivity1_Ext cOrder1_Con cDiscipline1_Con if TrumpIndex>.3920479 & TrumpIndex~=.

  set scheme stcolor
  estat ic
 
margins if cIdeas1_Open>.2594228  & cIdeas1_Open<.7833616 & cAesth1_Open>.2609113 & cAesth1_Open<.8215828  & cDepress1_Neu > .2305085  & cDepress1_Neu < .8531073 & cAnxiety1_Neu >  .0625  & cAnxiety1_Neu< .875 & cComply1_Agr > .3320099  & cComply1_Agr < .9434242   &  cAltruism1_Agr>.3446154 & cAltruism1_Agr<.8853846   & cAssertive1_Ext>.2 & cAssertive1_Ext<.8 & cActivity1_Ext>.2580645  & cActivity1_Ext<.8670968  & cOrder1_Con>.3736842  & cOrder1_Con < .8947368 & cDiscipline1_Con>.3313725 & cDiscipline1_Con < .9624183 , dydx( *) 

  marginsplot,  title(Trump Status) subtitle("") ytitle("") derivlabels xtitle("") horizontal plotopts(connect(none))  xline(0) xscale(r(-.6 .6)) xlabel(-.6(.2).6) xtick(-.6(.2).6)
  
     
graph save "Graph" "statusMean.gph", replace
 
  *****combine
  

graph combine jobMean.gph BigLieMean.gph leaderMean.gph statusMean.gph , rows(2) cols(2)  title(Aspects of Trump Followers) name(AspectsMean, replace) saving(AspectsMean, asis replace)
 
 
spearman Trump_JobApproval Trump_BigLie Trump_Leader Trump_Status if TrumpIndex>.3920479 & TrumpIndex~=., pw stats(rho obs p) star(0.01)
 
 

